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ABSTRACT 

Relativistic X-ray disk-lines have been found in multiple neutron star low-mass X-ray binaries, in close 
analogy with black holes across the mass-scale. These lines have tremendous diagnostic power and have been 
used to constrain stellar radii and magnetic fields, often finding values that are consistent with independent 
timing techniques. Here, we compare CCD-based data from Suzaku with Fe K line profiles from archival data 
taken with gas-based spectrometers. In general, we find good consistency between the gas-based line profiles 
from EXOSAT, BeppoSAX and RXTE and the CCD data from Suzaku, demonstrating that the broad profiles 
seen are intrinsic to the line and not broad due to instrumental issues. However, we do find that when fitting 
with a Gaussian line profile, the width of the Gaussian can depend on the continuum model in instruments 
with low spectral resolution, though when the different models fit equally well the line widths generally agree. 
We also demonstrate that three BeppoSAX observations show evidence for asymmetric lines, with a relativistic 
disk-line model providing a significantly better fit than a Gaussian. We test this by using the posterior predictive 
p-value method, and bootstrapping of the spectra to show that such deviations from a Gaussian are unlikely to 
be observed by chance. 

Subject headings: accretion, accretion discs — stars: neutron — X-rays: binaries 
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1. INTRODUCTION 

Broad iron emission lines are known to be prominent in the 
X-ray spectra of many accreting objects, especially in stellar- 
mass and supermassive black holes (for a review, see Miller 
2007). The highest quality X-ray observations suggest that a 
clear asymmetry is seen in the line profile, as would be ex- 
pected if the lines originate from the innermost region of the 
accretion disk and therefore subject to strong relativistic ef- 
fects. An origin in the innermost, relativistic regions of the ac- 
cretion disk is the favored explanation for the broad lines seen 
in black hole X-ray binaries and AGN (for example, see re- 
views bv lFabian et alj2000tlRevnolds & NowaU2003t [Miller 
120071 and references therein). This interpretation has held 
up to many alternative models for line bro adening and the na 
ture of the innermost accretion flow (e.g., Reynolds & Wilms 



neutron star LMXBs (Bhattacharwa & Strohmayer 2007; 



2010). 



2000; ReynoldsetaD 120091: iFabian et alJl2009t IMiller etal 



In neutron star low-mass X-ray binaries (LMXBs), iron 
emission lines were first discovered in t he mid-1980s with 
EXOS AT and Tenma (IWhite etal] 1 1981 \l9M. iHirano et all 
119871) and have been observed by every other major X- 
ray mission since (e.g. RXTE, ASCA, BeppoSAX, Chandra, 
XMM-Newton and Suzaku). Compared to the lines in black 
hole X-ray binaries and AGN, the lines seen in neutron star 
LMXBs are typically much weaker, with a peak deviation 
from the continuum of less than 10%. This likely hindered 
the detection of relativistic lines in neutron stars, which are 
otherwise expected since the potential around a neutron star 
is somewhat similar to that around a black hole with a low 
spin parameter. 

Recently, sensitive CCD-based spectroscopy with Suzaku 
and XMM-Newton has started to probe the line shape in 

ecackett@wayne.edu 



Cackett et al. 2008, 2009, 2010; Papitto et al. 2009; Reis et al. 
20091: ID' Ai et al JI20091 120 1 Ob iDi Salvo et dl2009tllaria et alj 
2009; Egr on et al.l 1201 ll) . In many cases, relativistic line 
shapes are found with inner disk radii, measured using rel- 
ativistic line models, close to expectations for the stellar ra- 
dius. Moreover, neutron star magnetic field strength esti- 
mates based on the inner disk radius for two acc reting pulsars 
are consistent with measurements from timing (|Cackett et alj 
l2lMlMliIer et al.11201 UlPapitto et al .1120091 1201 11) . Gas spec- 
trometer obs ervations have also revealed relativ istic lines in 
4U 1705-44 dPiraino et al.ll200HlLin et al.ll2010h . 

For bright sources, phot on pile-up can affect CCD spectra 
dBalletlll999l iDavisl feOQl). Pile -up occurs when more than 
one photon arrives in an event box in a single detector frame 
time. Thus, rather than recording the energy and time of each 
individual photon, all photons arriving at that location within 
a single frame time get counted as a single photon, with an 
energy e qual to the sum of the energy o f the individual pho- 
tons (see lDavisll200 It IMiller et alj |2010, for a more detailed 
discussion of this effect). Therefore, the shape of the detected 
spectrum is altered, with the spectrum becoming artificially 
stronger at higher energies and weaker at lower energies. This 
could, of course, alter the ability to accurately recover the iron 
K emission line profile in observations s trongly affected by 
this process. In order to address this issue. [Miller et al.l (1201 Oh 
performed a wide range of simulations based on typical black 
hole and neutron star LMXB spectra. These authors found 
that pile-up does not broaden the line profiles recovered, and 
in the most extreme cases it acts in the opposite sense, leading 
to a narrower line profile. 

An anal ysis of XMM-N ewton spectra of neutron star 
LMXBs by iNg et al.l (120101) came to an opposite conclusion 
about how pile-up affects the line profiles. These authors 
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looked at XMM-Newton spectra taken in pn timing mode, 
where the detector is read out continuously, and found that 
pile-up was significant in many of the observations. To miti- 
gate these effects, they extracted spectra only from the wings 
of the point-spread function (PSF), excluding the columns on 
the detector where the count rate was highest and thus has the 
most severe problems. They find that the resulting spectra are 
consistent with a Gaussian line profile, and that relativistic, 
asymmetric line models are not required to fit the data. These 
results differ from those obtained with Suzaku, wherein excis- 
ing successive annular regions of t he PSF yields highly con- 
sistent and relativistic l ine profiles (Cackett et al. 20 10J). Al- 
though iNg^eta]] (120 lOt) note that continuum modeling, back- 
ground subtraction, charge transfer inefficiencies, and X-ray 
loading of offset maps (an instrumental effect that is unique 
to the EPIC-pn came ra) have a significant effect on the over- 
all line profile, Ng et al. (2010) conclude that photon pile-up 
falsely leads to relativistic line profiles. However, in the spec- 
tra excluding the core of the PSF, the quality of the spectra 
are severely limited as the vast majority of source counts are 
excluded. 

In order to place new results from CCD spectroscopy of 
Fe K lines in neutron stars into a broader context, and in or- 
der to better understand the results of data from CCD and 
gas spectrometers, we systematically compare data from three 
gas spectrometer missions that are unaffected by pile-up (EX- 
OSAT, RXTE, BeppoSAX) with CCD data from Suzaku. We 
choose to compare with the Suzaku d ata, because, as dis- 
cussed in length by Mill er et al.l ((2010), pile-up is much less 
of a problem with the Suzaku/XIS detectors than with XMM- 
Newton/pn - the much broader PSF of Suzaku means that the 
source flux is sprea d over many more pixels. As shown by 
ICackett et al.l ((2010), extracting the source spectrum from an 
annulus, excluding the central core of the Suzaku PSF, does 
not change the shape of the line profile obtained. In this paper, 
we show that the gas spectrometer data give good consistency 
between the different missions and with the CCD Suzaku data, 
supporting the relativistic interpretation of the iron line origin. 

2. DATA SAMPLE AND REDUCTION 

In this paper, we analyze data from four different missions: 
EXOSAT, RXTE, BeppoSAX, and Suzaku. The detectors we 
utilize from EXOSAT, RXTE and BeppoSAX are gas spectrom- 
eters, whereas the Suzaku detectors covering the Fe K line re- 
gion are CCDs. The combination of both effective area and 
energy resolution is important for iron line spectroscopy and 
for the discussion of our results. The parameters for the detec- 
tors we use follows here. The EXOSAT/GSPC has an energy 
reso lution (FWHM) of approximately 10% at 6 keV (-0.6 
keV lPeacock et al.1 fl981) and an effective area of around 150 
cm 2 at 6 keV. The BeppoSAX/MECS instrument is very com- 
parable to EXOSAT/GSPC, with a simil ar energy resolutio n 
(-0.5 keV) and effective area at 6 keV dBoella et al.lll997l) . 
The proportional counters on RXTE, however, have a sig- 
nificantly higher effective area (about 1300 cm 2 per PCU), 
but a significantly lower energy resolution o f approximately 
FWHM ~1 keV at 6 keV (iJahoda et al.l 19961) . By comparison 
the two working front-illuminated XIS detectors on Suzaku 
have a combined effective area of around 600 cm 2 at 6 keV, 
and by far the best energy resolution of the detectors cons id- 
ered here, at approximately 0. 13 keV (Mit suda et al.ll2007b . 

Broadband energy coverage is also important in determin- 
ing the shape of the continuum either side of the iron line 



region. The energy range of EXOSAT/GSPC is dependent on 
the gain used during the observation. Most of the observations 
considered here are restricted to the 2-16 ke V range (Gain 1 ) , 
though some extend to higher energies (Gain 2), but in almost 
all those cases the S/N drops sharply above 20 keV. For RXTE 
we analyze data from the PCA only, which typically gives 
good spectra in the range from 3-30 keV. BeppoSAXhas sev- 
eral detectors to provide a broad energy coverage. Here, we 
used the LECS (0.5 - 3.5 keV), MECS (2 - 10 keV) and PDS 
(15 - 200 keV, with the upper limit dependent on S/N) to pro- 
vide broad energy coverage. In the case of Suzaku, we fit the 
XIS detectors in the energy range 1-10 keV, and use the PIN 
hard X-ray detector to provide broad band coverage from 12 
- 50 keV (with the upper limit dependent on S/N). 

We restrict our choice of sources here to four sources that 
have consistently shown the strongest broad iron emission 
lines in multiple observations by different missions, namely 
Serpens X-l, GX 349+2, GX 17+2 and 4U 1705-44. We 
analyze every EXOSAT and BeppoSAX observations of these 
sources. For Suzaku we use the Suzaku spectra analyzed by 
ICackett et al.l (|2010) here, please see that paper for details of 
the data reduction of those data. With RXTE many observa- 
tions of each source exist. As the focus of this work is on 
a comparison of the line profiles from different missions, as 
opposed to a detailed study of the line variability in each of 
the sources, we choose to look at only 3 RXTE spectra of each 
source. The specific spectra were initially chosen to be the 
three single longest continuous spectra of each source, but to 
ensure that we sampled a range of source states, we chose 
other observations in some cases. See Table Q]for details of 
the observations analyzed here. 

For EXOSAT, and BeppoSAX we used the pipeline pro- 
duced spectra available from the HEASARC database. The 
EXOSAT/GSPC pipeline produced spectra are time-averaged, 
background subtracted spectra. For bursting sources, the 
spectra have bursts already removed. The pipeline produced 
spectra are provided along with the corresponding response 
matrices. Further details of the EXOSAT data archive are 
given on the HEASARC websit^B The BeppoSAX pipeline 
spectra for the LECS, MECS and PDS are used. The PDS 
data are already background-subtracted, for the LECS and 
MECS we use the appropriate background file for the 8 ar- 
cmin extraction radius used for the source spectra, along 
with the appropriate response matrices from Se ptember 1997. 
For Su zaku, see the extensive details given in Cack ettet al.l 
(I20TOI) . here we use the spectra from that paper directly. Fi- 
nally, for RXTE, we use spectra from PCU2 only to provide 
the most reliable spectrum, with the data reduced following 
the standard procedures. Briefly, the standard 'goodtime' and 
'deadtime' corrections are applied, selecting data when the 
pointing offset < 0.02° and the Earth-limb elevation angle 
was > 10°. Spectra from PCU2 were extracted from the Stan- 
dard 2 mode data, and 0.6% systematic errors were applied to 
every channel. The appropriate background model for bright 
sources was used and response matrices were created with the 
PCARSP tool. 

3. SPECTRAL FITTING AND ANALYSIS 

Spe ctral fitting is performed with XSPEC vl2 (lArnaudl 
1996), and uncertainties are quoted at the la level of confi- 
dence throughout. 

1 http://heasarc.gsfc.nasa.gov/docs/joumal/exosat_archive6.html 
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TABLE 1 

Details of observations analyzed 



Source 


Mission 


Obs. No. 


Obs. ID 


Start date 
(dd/mm/yy) 


Exp. time 
(ksec) 


Source state 


References 


Serpens X.- 1 


FY/~IV A T 
CA fJ/i l 




6A9A« 


AX/AO/R 1 ; 

Uo/uy/ttj 


50 5 


Han ana 


12 3 4 


Serpens X- 1 






ZUo_i JUUZ 


AS/AO/QQ 

uj/uy/yy 


1 S 7/^9 0/1 S A 
IJ.1 1 dLXM 1 J.4 


Banana 




Serpens X-l 


RXTE 


1 


0fVV70 ai ni Ann 

ZUU / Z-U 1 -U 1 -UUU 


1 Q /A7 /Q7 

is/u //y / 


1 Z.o 


Banana 


0, / 


Serpens X-l 


RXTE 


Z 


iivm m O/i aaa 

ZUU / Z-U 1 -U4-UUU 


-i 1 ten /qi 
j i/u //y / 


14.3 


Banana 


0, / 


Serpens X- 1 


RXTE 


3 


40426-01-01-000 


05/09/99 


15.4 


Banana 


6.7 


Serpens X. - 1 


Sii ^xikit 


1 


401048010 


94/1 0/06 


18/29 


Banana 


8.9 


UA _14y+Z 


F YD C /I T 


I 


34067 


i a/ao/wj. 
i u/uy/o4 


18 8 


NB/FB 


2 10 


UA ,54y+Z 


CA L/J/i J 


■~) 


59605 


3 1 /AS /S 

3 1/Uo/oJ 


65 


NB/FB 


2 10 11 


ua j4y+z 


JJYD VA T 
CA J/1 i 


^ 


59663 


A1 /AO /s 1 ; 
u i/uy/oz> 


94 3 


MR /PR 


2 10 11 


ua j4y+z 


DCppOoft-A 




9 i oaqao i 
z i uuy uu i 


1 O/A'i/OO 
1 U/Uj/UU 


1 S IAA Q /9 9 S 

i j.y/44.y/zz. j 


PR iifitli MR/PR wf»i-i-,=»Y 

rrs wun iNo/ris venex 


12 


r.V -340-1-9 


Dcppc>3/\j\ 


2 


9 1 94.0009 

Z 1 ZH-UUUZ 


1 9/09/01 

1 Z/ UZ/ U 1 


1 5 6/76 1 /30 1 


FR wi\\\\ NFR/FR \ff=>i-fi=>Y 

FD W1L1I 1 > D/lD VCllCX 


1 3 


lja j4y+z 


DeppOolKA. 




9 1 9/10009 1 
Z 1 Z4UUUZ 1 


1 7 /A9 /0 1 
1 / /UZ/U 1 


99 ll'H'X AIA\ 
ZZ. //oJ.4/4 1 .U 


FR nfith MR/FR wf»i-i-,=»Y 

rrs wun iNo/rc venex 


i 


ua Jty+z 


RXTE 




3004.9 09 01 000 


00/01 /QS 


17 


FB 


14 


f.v -340-1-9 

UA JfH-y+Z 


RXTE 




30049 09 01 04 


1 3/01 /OS 


1 3 6 


FB 


14 


ua jn-y+z 


RXTE 




30049 09 09 060 


94/01 /OS 


1 5 7 


1MR/FR iipiTpy 


14 


ua J4y-t-z 


Suzaku 


1 


40000^010 
4UUUUJU 1 U 


1 zL/A'i/AA 
14/U_VUO 


s; /->() 


NB 


8 9 


UA _14y+Z 


Suzaku 


z 


4UUUUJUZU 


T O/AI /A6 

iy/uj/uo 


Q/Z4 


NB 


y o 


GX 17+2 


irv/~i cat 1 
CA UjA 1 


1 


337 1 5 


uj/uy/o4 


25. 1 


HB 


2, 15 


GX 17+2 


FY(1 VAT 


2 


33781 


06/09/84 


26. 1 


UR/IMR vprtPY 

I 1 Lj/ INJj VCllCA 


2 15 


GX 17+2 


CA (yj/i J 


^ 


58809 


90/08/8^ 
ZU/Uo/oJ 


6« 


NB 


2 15 


GX 17+2 


C A CM 1 


4 


60698 


i ^/ao/s 1 ; 
i z>/uy/oz> 


42 6 


HB 


2 15 


GX 17+2 


Ttorni/i^ A Y 

oeppOo/iA 




9 1 0^7001 
Z 1 U J / UU 1 


04/1 O/QQ 

U4/ lu/yy 


I 1 6//11 1/10 6 

I I .u/4 1.1/ iy.o 




1 6 


GX 17+2 




2 


9 1 0S7A01 1 
Z1UJ /UU1 1 


OS/1 O/QQ 

uj/ lu/yy 


99 0/7S 9/^6 Q 

zz.u/ / j.z/jo.y 


HR/MR 


1 6 


GX 17+2 


D ,,,,,,,, C/l V 


"2 

J 


9 1 O^TOOI 9 
Z1UJ /UU1Z 


07/1 O/QQ 

u // lu/yy 


1 J.Z/40. J/ZJ.4 


FTR/MR 


1 6 


GX 17+2 


Rfnnn^A Y 


4 


90961 01 1 

ZUZUi Ul 1 


03/04/07 


9 3/10 6/4. S 
L.Dl 1U.O/4.J 


HB 


| 7 


GX 17+2 


o€ppOo/iA 




90961 01 9 
ZUZ01 U1Z 


9 1 /A4/Q7 

z i/u4/y / 


1 476 4/9 6 
1 .4/O.4/Z.0 


HB 


| n 


GX 17+2 


RXTE 




20053-03-02-010 


07/02/97 


13.7 


NB 


18 


GX 17+2 


RXTE 


2 


30040-03-02-010 


18/11/98 


13.6 


FB 


18 


GX 17+2 


RXTE 


3 


30040-03-02-011 


18/11/98 


16.2 


FB 


18 


GX 17+2 


Suzaku 


1 


402050010 


19/09/07 


5/15 


NB 


9 


GX 17+2 


Suzaku 


2 


402050020 


27/09/07 


6/18 


NB 


9 


4U 1705-44 


BeppoSAX 


1 


21292001 


20/08/00 


20.6/43.5/20.1 


Soft 


19, 20 


4U 1705-44 


BeppoSAX 


2 


21292002 


03/10/00 


15.6/47.5/20.1 


Hard 


20,21 


4U 1705-44 


RXTE 


1 


20073-04-01-00 


01/04/97 


12.8 


Hard 




4U 1705-44 


RXTE 


2 


70038-04-01-01 


31/08/02 


12.6 


Soft 




4U 1705-44 


RXTE 


3 


90170-01-01-000 


26/10/05 


16.1 


Soft 




4U 1705-44 


Suzaku 


1 


401046010 


29/08/06 


14/14 


Hard 


9,21,22 


4U 1705-44 


Suzaku 


2 


401046020 


18/09/06 


17/15 


Soft 


9,21,22 


4U 1705-44 


Suzaku 


3 


401046030 


06/10/06 


18/17 


Soft 


9,21,22 



References. — (1) White etal. 1986, (2) Gottwald et al. 1995. (3) Schulzetal. 1989, (4) Seon & Min 2002, 
(5) Oosterbroek et al. 2001. (6) Muno et al. 2002. (7) Gladstone et al. 2007. (8) Cackettetal. 2008, (9) Cackett et al. 
2010, (10) Kuulkers & van der Klis 1995, (11) Ponman etal. 1988. (12) Di Salvo et al. 2001. (13) Iaria et al. 2004, (14) 
Izhang et al. 1998, (15) Kuulkers et al. 1997, (16) Di Sa lvo et all200ll (1 7)IFarinelli et aU2 005. (18) H oman et alf2002l 
(19) IPiraino et ali2007l (20) IFiocchi et alJ2007l (21) ILin et alJ20ig (22) IReis et 81120091 

NOTE. — Obs. No. is the observation number that we use in order to identify the observations. The BeppoSAX exposure 
times are given for the LECS, MECS and PDS (in that order), while the Suzaku exposure times are the good time exposure 
for the individual XIS detectors and the HXD/PIN good time (in that order). Source state abbreviations: NB - normal 
branch, HB - horizontal branch, FB - flaring branch. 



Broadband spectral fitting of neutron star LMXBs is often 
degenerate - a combinati on of several different models usu- 
ally fit the data well (e.g., Barret 20 011). This is a well known 
problem t hat led to the Eas tern dMitsuda et al.lfT989h versus 
Western (IWhite et al. 1988j) model debate. Here, we choose 
to explore two widely used continuum models for neutron 
star LMXBs in order to assess the importance of continuum 
choice on the iron emission line profiles. We first use a model 
consisting of a single-temperature blackbody, disk blackbody 
and power-law (i f needed), which we have successfully used 
in the past (e.g., ICackett et all 120081 120101) and is based on 
the re sults of testing multiple continuum models by dLin et al.l 
120071) . We refer to this model as Model 1. The second model 
we test consists of a single-temperature blackbody, Comp- 
tonized component (comptt in xspec), and in a number of 
cases an additional power-law (if needed). Hereafter, we refer 
to this model as Model 2 . This model i s widely used in the lit- 
erature for instance, see iBarretl (12001b ; iDi Salvo et aD (120001 
l20M : ltoiTetan(l200l . 

In both models Galactic photoelectric absorption is in- 



cluded with the phabs model. We mostly leave this as a 
free parameter in the model. However, with the RXTE data 
we alw ays fix the Nn parameter at the iDickey & Loc kman 
(1990) values (using the HEASARC N H tool), except for 
GX 17+2 where we use the A^h value determined from fitting 
the neutral absorption ed ges present Chandra gratings data 
dWroblewski et al.l 12008b, whic h is found to be higher than 
the Dickev & Lockman (1990) value. For spectra from the 
other missions, only when the A^h parameters tends to 0, or 
becomes significantly higher than the values in the literature 
do we fix the TVh- 

Initially, we fit the iron line with a Gaussian. All parameters 
(line energy, width and normalization) are free parameters in 
the fits, though the line energy of the Gaussian is restricted to 
be within the 6.4 - 6.97 keV range. 

Brief notes on individual observations follow. Please see 
Cackett et al. (2010) for notes on the Suzaku data. 

3.1. EXOSAT 

Ser X-l: The spectrum is fit from 2-16 keV. For 
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the fits with Model 2 we find an jVh significantly higher 
than the Dickey & Lockman (1990) value, thus for Model 
2 we fix A^h = 4 x 10 21 cm" 2 . This is similar to the 
value found when fitting XMM-Newton spectra of Ser X-l 
(Bhattach aryya & Strohmayedl2007l ICackett et aTll2010l) . 

GX 349+2: For all three observations we fit the data from 
2 - 20 keV (above 20 keV the S/N drops quickly). 

GX 17+2: We find for all observations that the spectrum 
turns up below 3 keV leading to too low an Nh (tending 
to zero). We therefore ignore below 3 keV and fix A^h = 
2.38 x 10 22 cm" 2 dWroblewski et al.ll2008h . this value is con- 
sistent with what we find when fitting the BeppoSAX and 
Suzaku observations of GX 17+2. The final observation (Ob- 
sID: 75122) has very poor calibration and there is a strong 
feature at 4.5 keV and between 10 - 14 keV (likely caused 
by variations in intrinsic detector gain, a known issue with 
EXOSAT). We therefore do not analyze that dataset here. For 
ObsID 33715 and 33781 we fit between 3-16 keV. For Ob- 
sID 58809 and 60698 a different gain was used, thus we fit 
from 3-20 keV (above 20 keV the S/N drops). 

4U 1705-44: We find no clear line detection in any of the 
observations. Upper limits on the line intensity are compa- 
rable or greater than the line intensities seen with other mis- 
sions, indicating that the spectra do not have the sensitivity 
to detect the Fe line. The strongest evidence for a line is in 
ObsID 60407 which is shown in Figure [TJ (bottom right hand 
panel, red). However, the Gaussian component in this obser- 
vation is required at less than the 2a level. We therefore do 
not consider these data any further here. 

3.2. RXTE 

Ser X-l: We fit the spectra from 3-20 keV. 

GX 349+2: The spectra are fit from 3-30 keV. 

GX 17+2: The spectra are again fit from 3-30 keV. 

4U 1705-44: The spectra are fit from 3-25 keV for obser- 
vation 2 and 3. For observation 1, where the source is in a 
hard state, we fit from 3-50 keV. 

3.3. BeppoSAX 

We fit the BeppoSAX spectra allowing an offset between the 
LECS, MECS and PDS spectra. The constant between LECS 
and MECS was constrained to be within 0.7 - 1.0, and the 
constant between the PDS and MECS was constrained to be 
within 0.7 - 0.95. Unless otherwise stated, we fit the LECS 
between 1-3.5 keV, the MECS from 2-10 keV and the 
PDS from 15 keV with the upper limit determined by the 
quality of the data at higher energies. As noted below, in 
several of the spectra, there is a feature at ^2.6 keV. This 
fe ature as been disc usse d in previous analy ses of these data 
bv llaria et all (120041) and lPiraino elaTI (120071) . who fit the fea- 
ture and discuss its potential origin. While the feature may be 
real, it occurs at a region in the spectrum of both the LECS 
and MECS where there is a significant change in the response 
of both detectors, and thus potentially instrumental. When 
we fit the feature with a Gaussian, we get paramet ers for this 
line th at are consistent wi t h the previous analysis of llaria et al.l 
j200l and iPiraino etall (I2007I) . Here, we choose to ignore 
this section of the data in the few cases where it is present, as 
it is not the focus of this work. This approach, or the approach 
of fitting it with a Gaussian lead to consistent spectral fits at 
better than the la level. We generally only fit the LECS data 
above 1 keV as we often found significantly reduced signal- 
to-noise there. 



Ser X-l: Here, the LECS data below 1 keV is of good qual- 
ity, thus we extend the spectral fitting range from 0.5 - 3.5 
keV. We fit the PDS from 15 - 30 keV. 

GX 349+2: We find a poor reduced chi-squared when fit- 
ting the standard MECS range due to a strong feature at 
around 2.6 keV in all 3 spectra. We therefore fit the MECS 
data from 3-10 keV in all cases. For the PDS data we fit the 
rang e 15-90 keV for ob servation 1 which has a strong hard 
tail (iDi Salvo et alj [2000). where as for observations 2 and 3 
where this hard tail is not significant (Iaria et al. 2004), we fit 
from 15-40 keV. We did not achieve any fits with x 2 < 2 
when using Model 1 . 

GX 1 7+2: We fit the PDS from 15-50 keV for all 5 obser- 
vations. 

4U 1705-44: In observation 1 we again find significant 
residuals at around 2.6 keV in the LECS and MECS spec- 
tra. We therefore fit the LECS data from 1 - 2.5 keV and the 
MECS data from 3-10 keV. We fit the PDS data from 15 - 
50 keV. Observation 2 displays a strong hard tail out to 100 
keV and thus we fit the PDS data up to 100 keV. However, an 
unbroken power-law does not fit the data well, thus we use an 
exponentially cut-off power-law for this observation. 

4. RESULTS 

Spectral fitting results are given in Tables |2]to|9] using con- 
tinuum models 1 and 2 and modeling the iron line with a 
Gaussian. We summarize the resulting line profiles in Fig- 
ure [TJ where we show a comparison of the iron line profiles 
for the four different objects and four different missions. For 
each object and mission we show only one spectrum. The 
profiles shown here are obtained from fitting the continuum 
model excluding the 5-8 keV range, and are plotted as the 
ratio of the data to the continuum model. We have chosen 
the continuum model that gives the best (lowest x 2 ) fit- The 
profiles show remarkable consistency between the different 
missions. All observations analyzed here give a positive de- 
tection of the iron, though this is no surprise given that these 
sources were chosen for that specific reason. We now discuss 
the main findings resulting from our analysis. 

4. 1 . Asymmetric line profiles 

Ou r previous analysis of Suzaku data (ICackett et al.ll2008l 
2010) has shown evidence for asymmetric line profiles, which 
can be fit well by a relativistic disk line model. The 
iron lines in EX OSAT and BeppoSAX have been analyzed 
in the past (e.g.. [Gottwald et alj|1995t IDi Salvo et all 120011: 
Oosterbroek et al. 2001; Iaria et al. " 120041) . however, gener- 
ally only a Gaussian line profile has been considered. An 
exception to this is the case of 4U 170 5-44, where both 
IPiraino et ail d2007l) and iLin et al.l (120101) have fitted a rela- 
tivistic diskline model to BeppoSAX data, finding that it is a 
better fit than a simple Gaussian. Here, we also fit the iron 
lines with a relativistic disk-line model (diskline, Fabi an et al.l 
1989). Of the eight EXOSAT spectra we analyze here, six 
show a decreased value of x 2 when using a diskline model 
rather than a Gaussian. However, in all those cases the fit 
with a Gaussian already provides a reduced-x 2 < 1, thus the 
improvement is not statistically significant. The parameters 
of the diskline fits are r easonable, and cons istent with typical 
parameters (such as in ICackett et al.l 1201 01) . but are not very 
well constrained. In the case of BeppoSAX, a diskline gener- 
ally provides equivalently good fits to the data. However, in 
three cases, we find that a diskline fits the data significantly 
better than a Gaussian. 
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FIG. 1. — A comparison of iron emission line profiles from the four different missions considered here, and good agreement is seen in most cases. Panel (a) 
shows Seipens X-l, (b) GX 349+2, (c) GX 17+2, and (d) 4U 1705^14. In all panels, black is SuzakufXIS (all FI combined), red is £XOSA77GSPG blue is 
BeppoSAXfMECS, and green is RXTE/PCA. 



TABLE 2 

Spectral fits from EXOSAT observations with the an absorbed diskbb + bbody + Gaussian model. 



Source 


Obs. No. N R Disk blackbody 


Blackbody 




Gaussian 


X 2 (dof) 




10 22 cm" 2 kT(keV) norm 


kT (keV) norm 


E (fceV) 


a (keV) nor 


n EW (eV) 



Serpens X-l 


1 


0.40 ± 0.15 


1.33 ±0.07 


61 ± 13 


2. 10 ±0.03 


(3.0±0.1) X 10" 


-2 


, ,.+0.07 
6 - 43 -0.03 


0.72 ± 0.09 


(6.9 ± 1.2) X 10" 


-3 


169 ± 36 


192.0(183) 


GX 349+2 


1 


0.25 ± 0.06 


2.09 ±0.10 


44 ± 7 


2.65 ±0.16 


(4.1 ±0.1) X 10" 


-2 


6.68 ± 0.05 


0.47 ± 0.08 


(1.4 ±0.2) X 10" 


-2 


102+* 6 


127.9 (127) 


GX 349+2 


2 


0.60 ± 0.04 


2.35 ± 0.05 


29 ± 2 


3.46 ± 0.27 


(2.5 ±0.5) X 10" 


-2 


640+ 0.08 


1 .09 ± 0.09 


(3.6 ±0.4) X 10" 


2 


221_ 23 


130.4(115) 


GX 349+2 


3 


0.75 ± 0.05 


2.44 ± 0.03 


29 ± 1 


4.39 ± 0.45 


(1.9 ±0.1) X 10" 


-2 


6.63 ±0.10 


0.85 ±0.15 


(2.5 ±0.4) X 10" 


2 


155 ± 16 


120.1 (115) 


GX 17+2 


1 


2.3S (fixed) 


1.40 ±0.03 


192 ± 16 


2.57 ± 0.02 


0.117±0.002 




6.58 ±0.10 


0.80 ±0.14 


(1.9 ± 0.5) X 10" 


2 


131 ±23 


139.6(185) 


GX 17+2 




2.3S (fixed) 


1.55 ±0.04 


151 ± 14 


2.46 ± 0.03 


(9.2 ±0.4) X 10" 


-2 


6.68 ±0.14 


0.80 ± 0.26 


(1.2 ±0.6) X 10" 


2 


81 ±20 


157.8 (186) 


GX 17+2 


3 


2.3S (fixed) 


1.67 ±0.03 


103 ±6 


2.54 ±0.04 


(5.0 ±0.2) X 10" 


-2 


6.74 ±0.08 


0.54 ±0.11 


(1.2 ±0.2) X 10" 


2 


99 ± 15 


97.7(112) 


GX 17+2 


4 


2.38 (fixed) 


1.46 ±0.03 


144 ± 11 


2.80 ±0.02 


0.1 16 ± 0.001 




64+ 0.16 


1.87 ±0.14 


(5.9 ± 0.9) X 10" 


2 


.,,+329 


202.0(124) 



TABLE 3 

Spectral fits from EXOSAT observations with the an absorbed bbody + comptt + Gaussian model. 



Source 


Obs. No. N H Blackbody 


Comptonized component 




Gauss 




X 2 (dof) 




10 22 cm" 2 kT(keV) norm 


T Q (keV) kT (keV) t no 


■m E (keV) 


cr (keV) 


norm 


EW (eV) 



Serpens X— 1 


1 


0.40 (fixed) 


nH7 +0.19 
°- 87 -0.05 


(1.1 ± 0.4) 




10 


-2 


15+°' 1 1 
UJ3 -0.35 


2.22 ± 0.06 


8.6+J- 6 


057 +0.14 
U3/ -0.04 


6.42+ 009 


0.73 ±0.10 


(7.0 ± 1.7) X 10' 


3 




191.6(182) 


GX 349+2 


1 


0.58 ± 0.32 


1.33 ±0.08 


(5.1 ± 0.4) 




10 


-2 


0.46 ± 0.05 


2.47 ± 0.06 


7.9 ± 0.5 


Sffl 

-0,1 


6.71 ±0.05 


0.39 ± 0.08 


(1.2±0.2)X 10" 


2 


124.6(125) 


GX 349+2 




Q.5(fixed) 


1 .44 ± 0.02 


(9. 1 ± 0.3) 




10 


-2 


33+ - 05 
U ---0.32 


2.91 ±0.06 


6.8 ± 0.2 


6.84 ± 0.05 


0.31 ±0.07 


(9.7 ± 1.4) X 10" 


3 


67 ± 9 


101.1 (113) 


GX 349+2 


3 


54+ - 55 
U ' 3 -0 .17 


1.52 ±0.03 


(9.9 ± 0.3) 




in 


-2 


0.47 ±0.10 


3.15 ±0.12 


6.0 ± 0.5 


9 +(>i 
Uy -0.2 


6.82 ± 0.06 


0.34 ±0.07 


(1.0 ±0.2) X 10" 


2 


64+1 7 


109.9(114) 


GX 17+2 


1 


2.38 (fixed) 


1.11 ±0.09 


(2.5 ± 0.6) X 10 


-2 


0.43 ± 0.09 


3.09 ± 0.08 


6.0 ± 0.3 


1.6 ±0.2 


6.67 ±0.12 


0.48 ±0.16 


(8.8 ± 3.1) X 10" 


3 


61 ± 16 


120.9 (183) 


GX 17+2 


2 


2.38 (fixed) 


1.32 ±0.08 


(5.0 ± 0.6) 




in 


-2 


0.51 ± 0.05 


2.96 ±0.14 


5.6 ± 0.4 


1.5 ±0.1 


6.75 ±0.12 


0.28 ± 0.23 


4.6^;| X 10" 3 




31 ± 12 


142.6(184) 


GX 17+2 


3 


2.38 (fixed) 


1.19 ±0.02 


(5.8 ± 0.2) 




in 


-2 


0.25 ± 0.06 


2.68 ± 0.04 


6.1 ± 0.1 


2.3 ±0.5 


6.79 ± 0.06 


0.44 ± 0.09 


(9.4 ± 1.7) X 10" 


3 


79 +175 


100.8(117) 


GX 17+2 


4 


2.38 (fixed) 


1.46 ±0.09 


(1.7 ± 0.2) 




in 


-2 


0.24 ± 0.05 


3.34 ± 0.05 


5.8 ± 0.1 


2.2 ± 0.3 


6.79 ± 0.09 


0.41 ±0.11 


(7.2 ± 1.5) X 10" 


3 


55+ 248 


150.4(129) 
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TABLE 4 

Spectral fits from RXTE observations with the an absorbed diskbb + bbody + power-law + Gaussian model. 



Source 


Obs. No. 




Disk blackbods 




Blackbody 


Power- law 


Gaussian 






X 2 (d°f) 






10 22 em- 2 


kT (keV) 




kT (keV) 




T norm E (keV) 


a (keV) 






EW (eV) 




Serpens X- 1 


1 


0.4 (fixed) 


1.54 ±0.03 


49 ± 3 


2.37 ± 0.02 


(2.3 ± 0.1) X 10" 2 


6.51 ±0.06 


0.48 ±0.12 


(4.8 ±0.7) X 10" 


-3 


104 ± 21 


27.0 (38) 


Serpens X- 1 


2 


0.4 (fixed) 


1.66 ±0.04 


50 ±4 


2.39 ± 0.03 


(3.3 ± 0.2) X 10" 2 


6.51 ±0.06 


0.64 ± 0.11 


(9.6 ± 1.5) X 10" 


-3 


138 ± 18 


25.1 (38) 


Serpens X- 1 


3 


0.4 (fixed) 


1.71 ±0.04 


41 ±3 


2.41 ± 0.03 


(2.9 ± 0.2) X 10" 2 


6.51 ±0.06 


0.52 ±0.12 


(7.0 ± 1.1) X 10" 


-3 


110 ±20 


18.1 (32) 


GX 349+2 


1 


0.5 (fixed) 


1.98 ±0.03 


69 ± 3 


2.64 ± 0.03 


(6.5 ± 0.3) X 10" 2 


6.40+ - 01 


1.05 ± 0.05 


(5.9 ±0.3) X 10" 


-2 


302 ± 23 


46.2 (51) 


GX 349+2 


2 


0.5 (fixed) 


2.28 ± 0.04 


46 ± 3 


2.87 ± 0.06 


(5.1 ±0.6) X 10 -2 


6.40+0-02 


1.26 ±0.06 


(6.8 ± 0.4) X 10" 




306 ±25 


53.4 (52) 


GX 349+2 


3 


0.5 (fixed) 


1.77 ±0.02 


87 ±4 


2.56 ± 0.02 


(6.4 ±0.2) X 10" 2 


6.40+0.07 


0.82 ±0.09 


(2.9 ±0.3) X 10" 


-2 


186 ± 18 


37.1 (51) 


GX 17+2 


1 


2.38 (fixed) 


1.70 ±0.01 


130 ±4 


2.68 ± 0.01 


(7.5 ± 0.2) X 10" 2 


6.51 ±0.10 


0.64 ±0.13 


(1.5 ±0.3) X 10" 


-2 


86 ± 14 


45.2 (51) 


GX 17+2 


2 


2.38 (fixed) 


2.05 ± 0.03 


81 ±3 


2.76 ± 0.02 


0.125 ±0.004 


ft ,,,+0.09 
6 - 46 -Q.Q6 


0.85 ±0.12 


(3.9 ±0.6) X 10" 


-2 


135 ± 16 


38.5 (52) 


GX 17+2 


3 


2.38 (fixed) 


1.88 ±0.02 


105 ±4 


2.64 ± 0.02 


(9.7 ±0.3) X 10" 2 


640 +0.05 


0.92 ± 0.07 


(5.4 ± 0.4) X 10" 


-2 


218 ± 19 


49.0 (52) 


4U 1705-44 


1 


0.8 (fixed) 


2.28 ± 0.05 


2.7 ±0.2 


4.68 ± 0.01 


(1.2 ±0.2) X 10" 2 


0.6±0.8 (2.1+9[j) X 10~ 3 6.40+ 001 


0.62 ± 0.07 


(3.0 ±0.2) X 10" 


-3 


197 ± 49 


65.2 (75) 


4U 1705-44 


2 


0.8 (fixed) 


1.99 ±0.04 


18 ± 1 


2.62 ± 0.05 


(1.62 ±0.01) X 10" 2 


64Q +0.02 


0.57 ± 0.09 


(6.7 ± 0.7) X 10" 


-3 


130 ± 16 


30.6 (42) 


4U 1705-44 


3 


0.8 (fixed) 


1.52 ±0.04 


19 ± 2 


2.27 ± 0.02 


(1.18 ±0.05) X 10" 2 


6.52 ± 0.07 


0.32 ±0.15 


(1.6 ±0.3) X 10" 


-3 


82 ± 26 


33.8 (42) 



TABLE 5 

Spectral fits from RXTE observations with the an absorbed bbody + comptt + Gaussian model. 



Source Obs. ID JVjj Blackbody Comploni zed component Gaussian x (dof) 







10 22 cm- 2 


kT (keV) 


norm 




7" (keV) 


kT (keV) 




norm 


E (keV) 


a (keV) 


norm 




EW (eV) 




Serpens X-l 
Serpens X-l 


1 

2 


0.4 (fixed) 
0.4 (fixed) 


1.14 ±0.07 
1.21 ±0.08 


(1.3 ±0.3) X 10" 
(2.1 ±0.4) X 10" 


-2 


041 +0.14 

n 44+o:t§ 

u ' -0.43 


2.58 ± 0.04 
2.51 ± 0.04 


6.0 ±0.3 
6.8 ±0.3 


72 +28.08 

(\ a-j.+^ie 


6 45+0.08 

6 45+8:85 
"■ 43 -0.05 


0.53 ± 0.15 
0.66 ±0.14 


(5.2 ± 1.4) X 10" 
(9.9 ±2.7) X 10" 


-2 
-3 


no+'f 

138+1* 


22.0 (36) 
21.6 (36) 


Seipens X- 1 


3 


0.4 (fixed) 


1.23 ±0.06 


(2.0 ±0.2) X 10' 


-2 


39 +0.25 


2 50 +0.1 1 
l - 3 -0.03 


6.8 ± 0.3 


n s?13.80 
0.82_o 17 


6.47 ± 0.07 


0.53 ± 0.17 


< 7 °-f:o> x 10 ~ 


3 


,°<4 


15.7(30) 


GX 349+2 


1 


0.5 (fixed) 


1 36 +0 - 10 
lib -QM 


(7.5 ± 2.3) X 10" 


-2 


°- 59 -0.58 


2.64 ± 0.05 


6.8 ±0.8 


■ .+13.8 

i> 


640+ 0.20 


nn+0.09 
uw -0.15 


<4 9 ^2 ' X 10 


2 


246^ 


36.8 (49) 


GX 349+2 


2 


0.5 (fixed) 


1.42 ±0.03 


0.11 ±0.01 




0.40 ±0.37 


2.67 ± 0.03 


7.6 ±0.4 


1J1 -0.06 


ft 5S+0-13 

+8:15 

°- -0.03 


0.90 ± 0.28 


(3.0 ± 1.3) X 10" 


-2 


139+lf 1 


32.4 (43) 


GX 349+2 


3 


0.5 (fixed) 


J - -0.04 


(5.6 ± 2.1) X 10" 


-2 


56+°- 14 
U3t) -0.55 


2.65 ± 0.06 


6.6 ± 0.7 


0.73 ± 0.20 


(2.5 ± 1.0) X 10" 


-2 


162 ±45 


34.3 (49) 


GX 17+2 


1 


2.38 (fixed) 


1.83 ±0.25 


(1.9 ± 0.8) X 10" 


-2 


0.73 ± 0.03 


3.1 ± 0.2 


4.5 ±0.5 


1.6 ±0.1 


6.40+ 004 


1.03 ± 0.09 


(4.1 ±0.7) X 10" 


-2 


228+ 74 


28.9 (49) 


GX 17+2 


2 


2.38 (fixed) 


1.91 ±0.18 


(6.1 ± 1.3) X 10" 


-2 


0.72 ± 0.04 


3.0 ± 0.1 


5.5 ± 0.6 


2.0 ±0.2 


6.40+ - 07 


1.10 ± 0.10 


(7.1 ± 1.3) X 10" 


-2 


248™ 
318+J 


34.3 (50) 


GX 17+2 


3 


2.38 (fixed) 


1.79 ±0.24 


(4.3 ±0.9) X 10" 




0.75 ± 0.04 


2.9 ± 0.1 


5.2 ±0.6 


2.0 ±0.1 


640+ 0.05 


1.06 ±0.08 


(7.7 ± 1.0) X 10" 


-2 


42.4 (50) 


4U 1705-44 


1 


0.8 (fixed) 


1.72 ±0.10 


(1.7 ± 0.6) X 10" 


-3 


0.80 ±0.03 


10.7 ±0.9 


2.9 ±0.2 (2.: 


i±0.2) X 10~ 2 


6 . 4 -rt).02 


0.69 ± 0.08 


(3.4 ± 0.4) X 10" 


-3 


229 ± 38 


66.1 (75) 


4U 1705-44 


2 


0.8 (fixed) 


1 87+0.78 
2H7 -0.07 


(9.6+^) X 10 _: 


? 


0.73 ± 0.02 


2.0+0.3 


7.6 ±0.3 


0.62+ - 01 


6.40+ 006 


0.67 ±0.10 


(8.0 ± 1.0) X 10" 


-3 


158 ± 30 


27.4 (40) 


4U 1705-44 


3 


0.8 (fixed) 


1.17 ±0.08 


(3.5 ±0.7) X 10" 


-3 


32 +0.14 


2.4 ± 0.1 


6.8 ±0.3 


35+§ : ^ 
u -"-0 09 


6 45+ -(> 7 


0.34 ±0.15 


(1.6 ±0.4) X 10" 


-3 


89 +5559 


23.0 (40) 



TABLE 6 

Spectral fits from BeppoSAX observations with the an absorbed diskbb + bbody + power-law + Gaussian model. 



Source 


Obs. 


Nn Disk blackbody 


Blackbody 




Powe 


r-law 










X 2 (dof) 






10 22 cm" 2 kT(keV) 


norm 


kT (keV) 


norm 




r 


norm 


E(keV) 


cr (keV) 


norm 




EW (eV) 




Serpens X- 1 


1 


0.62±0.01 1.46±0.02 


51 ±3 


2.21 ± 0.04 


(3.0±0.1) x 10" 


-2 


2.98 ± 0.06 


0.90 ± 0.05 


6 40+0-01 


0.96 ± 0.08 


(1.2 ±0.2) X 10" 


-2 


234 ± 31 


222.2 (133) 


GX 349+2 


1 


No fit with reduced x^ < 2 


























GX 349+2 


2 


No fit with reduced % 2 < 2 


























GX 349+2 


3 


No fit with reduced x 2 < 2 












n 97+ y ---' 
u.z/_ . 19 














GX 17+2 


1 


2.12 ±0.04 1.81 ±0.03 


77 ± 5 


2.79 ± 0.02 


(7.0 ±0.3) X 10" 


-2 


2.4 ±0.2 


6.60 ± 0.06 


0.53 ±0.15 


(9.1 ± 2.2) X 10" 


-3 


60 ± 8 


203.1 (119) 


GX 17+2 


2 


2.12 ±0.05 1.87 ±0.01 


70 ± 2 


2.79 ± 0.02 


(5.9 ±0.2) X 10" 


-2 


2.7 ±0.3 


6.70 ± 0.03 


0.33 ±0.05 


(5.8 ±0.6) X 10" 


-3 


40 ±4 


231.5 (119) 


GX 17+2 


3 


No fit with reduced % 2 < 2 
























GX 17+2 


4 


2.27 ±0.10 1.70 ±0.06 


74 ± 11 


3.01 ± 0.02 


0.104 ± 0.005 




2.6 ±0.1 


1.6 ±0.6 


6.83 ±0.12 


0.97 ± 0.20 


(2.0 ±0.6) X 10" 




154 ± 48 


175.7 (120) 


GX 17+2 


5 


2.08±0.1O 1.72±0.06 


86 ± 12 


2.88 ± 0.04 


(9.8 ±0.6) X 10" 


-2 


2.3 ±0.6 


0.3+08 


6 50+ 016 
+0°28° 


0.72 ± 0.26 


(1.3 ± 0.6) X 10" 




84 ± 29 


143.7(120) 


4U 1705-44 


1 




42 ±3 




(2.9±0.1) X 10" 




2.8 ± 0.1 










191 ± 27 


162.6 (91) 


4U 1705-44 




1.48 ±0.04 2.03 ±0.07 


2.3 ±0.3 


4.4 ±0.2 


(5 ± 1) X 10~ 3 




c<- 2 * ,1 


■'S xl0 " 2 


6.60 ±0.11 


0.76 ± 0.24 


(1.4 ± 0.5) X 10" 




145 ± 88 


116.2 (106) 



Note. — * The second observation of 4U 1705-44 has a strong hard tail extending to 100 keV. A single power-law does not fit this tail well, thus we use a cut-off power-law instead. We find the cut-off energy, 
Ecut - 28 ± 5 keV. 



TABLE 7 

Spectral fits from BeppoSAX observations with the an absorbed bbody + comptt + power-law + Gaussian model. 









kT (keV) 






7- <keV) 


kT (keV) 






r 




E (keV) 








EW (eV) 




Serpens X-l 


1 


0.56 ±0.05 


0.99 ± 0.02 


(1.4 ± .1) x 10" 


-2 


0. 10 ± 0.01 


2.58 ± 0.04 


6.7 ±0.1 


2.0 ±0.2 






6.40+0-26 


1.01 ±0.07 


(1.5 ± 0.2) X 10" 


-2 


288+ 40 ) 


226.6(132) 


GX 349+2 


1 


0.52 ±0.03 


1.44 ±0.02 


(6.2 ± 0.3) X 10" 


-2 


0.45 ±0.01 


2.63 ± 0.04 


6.7 ±0.2 


1.18 =b 0.03 


2.1 ±0.3 


{4.S+W 2 ) X 10" 2 


6.68 ±0.03 


0.40 ± 0.04 


(1.1 ±0.1) X 10" 


-2 


80 ± 7 


160.7 (105) 


GX 349+2 


2 


0.52 ±0.01 


0.60 ± 0.01 


(6.2 ±0.1) X 10" 


-2 


1.39 ± 0.01 


2.94 ± 0.02 


4.2 ±0.1 


1.20 ±0.01 






6.76 ± 0.02 


0.20 ±0.03 


(5.9 ±0.4) X 10" 


-3 


31 ± 3 


182.9(100) 


GX 349+2 


3 


0.52 ±0.01 


0.62 ±0.01 


(6.0 ±0.1) X 10" 




1.43 ± 0.01 


3.08 ± 0.04 


3.6 ±0.1 


1.00 ± 0.02 






6.71 ±0.02 


0.24 ± 0.03 


(5.9 ±0.4) X 10" 


-3 


34 ± 3 


189.2(116) 


GX 17+2 


1 


2.27 ± 0.09 


1.50 ±0.05 


(6.6 ± 0.4) X 10" 


-2 


0.57 ± 0.01 


3. 10 ±0.07 


5.5 ±0.3 


0.99 ± 0.03 


3.1 ±0.1 


2.6 ±0.7 


6.68 ± 0.04 


0.21 ±0.05 


(4.0 ±0.6) X 10" 


-3 


27 ± 6 


146.7(117) 


GX 17+2 


2 


2.33 ±0.10 


1.44 ±0.03 


(6.8 ± 0.2) X 10" 


-2 


0.54 ±0.01 


3.01 ±0.06 


5.6 ±0.2 


1.06 ±0.02 


3.5 ±0.2 


3.0 ±0.8 


6.73 ± 0.03 


0.17 ±0.05 


(3.3 ± 0.4) X 10" 


-3 


23 ±4 


162.6(117) 


GX 17+2 


3 


2.18 ± 0.10 


1.54 ± 0.03 


(6.6 ±0.2) X 10" 


-2 


0.56 ± 0.01 


3.18 ±0.05 


5.5 ±0.1 


1.02 ±0.02 


3.2 ±0.1 


2.0 ± 0.6 


6.73 ± 0.05 


0.25 ± 0.08 


(3.1 ±0.6) X 10" 


-3 


21 ± 6 


174.7(118) 


GX 17+2 


4 


2.25 ±0.17 


1.59 ± 0.08 


(3.2 ±0.2) X 10" 


-2 


0.56 ± 0.03 


3.37 ± 0.06 


6.0 ±0.2 


0.89 ± 0.05 


2.8 ±0.1 


2.6± 1.1 


6.77 ± 0.06 


0.20 ±0.10 


(4.3 ± 1.0) X 10" 


-3 


32 ± 11 


150.1 (118) 


GX 17+2 


5 


2.10 ±0.29 


1.57 ±0.11 


(5.3 ± 0.7) X 10" 




0.57 ± 0.03 


3.25 ±0.14 


5.6 ±0.4 


0.99 ± 0.06 


3.0 ±0.3 


1.8 ± 1.5 


6.68 ±0.10 


0.17 ±0.17 


(3.3 ± 1.5) X 10" 


-3 


22 ± 15 


132.2(118) 


4U 1705-44 


1 


1.57 ±0.15 


1.51 ± 0.03 


(2.1 ± 0.2) X 10" 


-2 


0.44 ± 0.02 


2.69 ± 0.06 


5.9 ± 0.2 


0.49 ± 0.03 


2.8 ±0.2 


0.49 ±0.33 


6.56 ± 0.04 


0.27 ± 0.06 


(2.9 ± 0.4) X 10" 


-3 


58 ± 14 


141.9 (89) 


4U 1705-44 


2 


091 +0.28 








0.73 ± 0.01 


6.65 ± 0.52 


4.0 ± 0.1 


(3.2 ± 0.2) X 10" 2 


-0 74 + '-8 * 


(1.0+ 27 Q 2 ) X 10" 3 


6.4+ 04 


1.12 ±0.09 


(3.5 ±0.5) X 10" 


-3 


17S +81 


115.2(106) 



Note. — * The strong hard tail extending to 100 keV requires a cut-off power-law. We find the cut-off energy. E cu i = 20+Jj keV. 
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TABLE 8 

Spectral fits from Suzaku observations with the an absorbed diskbb + bbody + power-law + Gaussian model. 



Source Obs. % Disk blackbody Blackbody Power-law Gaussian x 2 (dof) 

10 22 cm" 2 kT(keV) norm kT(keV) norm r norm E (keV) ct (keV) norm EW (eV) 



Serpens X-l 


1 


0.60 ± 0.01 


1.24 ±0.01 


91 ±3 


2.17 ±0.02 


(5.1 ±0.1) X 10" 


2 


2.7 ±0.1 


0.65 ± 0.04 


6.65 ± 0.01 


0.26 ± 0.02 


(2.5 ±0.2) X 10 


-3 


44±3 


1567.7(1107) 


GX 349+2 


1 


0.79 ± 0.02 


1.60 ±0.03 


95 ±6 


2.25 ± 0.02 


0.117±0.003 




2.1 ±0.2 


0.17±0.11 


6.55 ± 0.03 


0.58 ± 0.04 


(1.4 d= 0.1) X 





-2 


87 ±6 


1341.9(1108) 


GX 349+2 


2 


0.88 ± 0.02 


1.50 ±0.02 


97 ±4 


2.30 ± 0.02 


(9.5 ± 0.1) X 10" 


-2 


2.4 ±0.1 


0.44 ± 0.09 


6.71 ±0.02 


0.35 ± 0.03 


(7.5 ±0.5) X 





-3 


62 ±4 


1293.6 (1 108) 


GX 17+2 


1 


2.24 ±0.11 


1.77 ±0.02 


83 ±4 


2.68 ±0.03 


(5.9 ±0.2) X 10" 


-2 


2.9 ±0.4 


0.84 ±0.68 


6.57 ±0.05 


0.28 ± 0.05 


(4.3 ± 0.7) X 





-3 


30 ±6 


644.0 (608) 


GX 17+2 


2 


2.27 ± 0.07 


1.88 ±0.02 


71 ± 3 


2.64 ± 0.02 


(7.6 ±0.3) X 10" 


-2 


3.0 ±0.2 


1.1 ±0.5 


6.60 ± 0.02 


0.21 ± 0.02 


(5.4 ±0.5) X 





-3 


32 ±4 


665.9 (608) 


4U 1705-44 


1 


2.0 (fixed) 


8.4 ±0.8) X 10" 2 


(7.7+l 7 g 2 ) X 10 7 


1.18 ±0.06 


(6.6 ±0.5) X 10" 


-4 


1.67 ±0.01 


0.127 ±0.02 


6.55 ±0.03 


0.20 ± 0.03 


(2.5 ±0.4) X 





-4 


41 ± 11 


961.8 (1123) 


4U 1705-44 


2 


2.01 ± 0.02 


1.23 ±0.02 


44 ± 3 


2.18 ±0.03 


<2.3±0.1)X 10" 


-2 


3.0 ±0.1 


0.79 ± 0.05 


6.4+ - 1 


0.68 ± 0.03 


(3.5 ±0.3) X 


10 


-3 


121 ± 16 


1364.5(1123) 


4U 1705-44 


3 


1.98 ±0.03 


0.74 ±0.01 


139 ± 15 


1.84 ±0.01 


(9.8 ±0.1) X 10" 


-3 


2.9 ±0.1 


0.38 ± 0.05 


6.4+ - 2 


0.40 ± 0.05 


(4.5 ± 0.8) X 


10 


-4 


40 ± 14 


1275.8(1122) 



TABLE 9 

Spectral fits from Suzaku observations with the an absorbed bbody + comptt + power-law + Gaussian model. 



Source 


Obs 


N H 




Blackbody 








Compto 


lized component 








G£ 








x 2 (dof) 






10 22 cm 2 


kT (keV) 








T (keV) 


kT (keV) 






r 


norm E (keV) 


o- (keV) 






EW (eV) 




Serpens X-l 


1 


0.57 ±0.01 


0.75 ±0.01 


(8.6 ± 0.4) 




10-3 


0.15 ±0.01 


2.44 ± 0.01 


7.8 ±0.1 


1.54 ±0.09 




6.65 ±0.01 


0.26 ± 0.02 


(2.5 ±0.1) X lO- 


-3 


44 ± 3 


1605.7(1109) 


GX 349+2 


1 


0.55 ±0.01 


1.37 ±0.07 


(4.8 ± 0.4) 




10" 2 


0.45 ± 0.01 


2.52 ± 0.04 


7.9 ±0.3 


1 .48 ± 0.04 




6.67 ± 0.03 


0.38 ±0.05 


IS. 6± 1.1) X 10" 


-3 


54 ±5 


1476.2(1108) 


GX 349+2 


2 


0.93 ±0.01 


0.94 ±0.01 


(3.00 ± 0.08 


) X 10" 2 


0.13 ±0.01 


2.51 ± 0.01 


8.31 ±0.06 


2.6 ± 0. 1 




6.65 ± 0.03 


0.49 ± 0.05 


(1.1 ±0.1) X 10 


-2 


59 ± 11 


1367.3(1108) 


GX 17+2 


1 


2.55 ± 0.06 


1.30 ±0.04 


(4.8 ± 0.4) 




KT 2 


0.60 ± 0.02 


2.83 ± 0.05 


6.1 ±0.2 


1.10 ±0.02 


3.5 ±0.1 


5.1 ±0.7 6.54 ±0.06 


0.31 ±0.07 


(4.6 ± 1.0) X 10 


-3 


32 ±7 


618.7 (607) 


GX 17+2 


2 


2.50 ±0.06 


1.36 ±0.05 


(5.3 ± 0.3) 




lO" 2 


0.56 ± 0.02 


2.73 ± 0.05 


6.8 ±0.3 


1.27 ±0.02 


3.4 ± 0.1 


4.6 ±0.6 6.60 ±0.02 


0.20 ±0.03 


(5.0 ±0.6) X 10 


-3 


30 ±4 


637.6 (606) 


4U 1705^14 


1 


1.56 ±0.03 


0.96 ± 0.01 


(2.23 ± 0.02) X 10" 3 


12+ - 12 
U1 -0.01 


6.6 ± 0.1 


5.94±0.03 (3.5 ± 0.1) X 10 -2 




64 +0.01 


0.87 ±0.05 


(1.4 ± 0.1) X 10 


-2 


-.ig+485 
Z -22 


1009.2(1123) 


4U 1705-44 


2 


1.90 ±0.01 


0.61 ±0.02 


(1.9 ±0.3) 




L<T 3 


0.14 ±0.02 


2.54 ± 0.01 


6.85 ± 0.07 


1 .0 ± 0. 1 




6.4+ - 2 


0.60 ± 0.03 


(2.7 ±0.2) X 10 


-3 


92 ± 15 


1368.2(1125) 


4U 1705-44 


3 


1.59 ±0.02 


1.98 ±0.03 


(7.1 ± 0.4) 




lO" 3 


0.36 ±0.01 


3.9 ± 0.5 


3.1 ±0.4 


0.13 ±0.02 




6.4+ 01 


0.46 ± 0.05 


(6.8 ± 1.0) X 10 -4 


62 ± 19 


1261.6(1123) 



The three cases where a diskline is a significantly better 
fit are GX 349+2 observation 1, GX 17+2 observation 1 and 
4U 1705-44 observation 1. Note that the 4U 1705-44 ob- 
servation is the same one discussed in detail by Pirai no et alj 
d2007l) . who also find an improvement in \ 2 when using a 
diskline compared to a Gaussian. The parameters for the 
diskline fits to these three spectra are given in Table [10] 
We give the parameters using the Model 2 continuum as 
these provide a better fit in all cases. A comparison be- 
tween the diskline fit parameters in Table [10] with previ- 
ously published fits to these sources show they are consis- 
tent. Our fits t o 4U 1705-44 are consistent with similar fits 
performed by iPiraino et al.l (|2007l) to the same data. Further- 
more, the diskline parameters we find here are all reasonable 
when comp aring with fits to Suzaku and XMM-Newton data of 
GX 3 49+2 dCackett et al.| l2008l I20T0I; llaria et al.ll2009]). GX 
17+2 dCackett et alj|2010f) and 4U 1 705-4 4 dReis et all 2009 1 : 
Pi Salvo et alj 120091: ICackett et all I20T0I iD'Ai et alj I2OI0I: 
Lin et alJ2010l) . In particular, the inner disk radius and equiv- 
alent width are consistent with this previous work for similar 
source states. 

In Figure [2] we show the line profiles with both a Gaus- 
sian and diskline model. The diskline fit to observation 1 of 
GX 349+2 gives an improvement of A\ 2 = 22.6 for 2 addi- 
tional degrees of freedom. This corresponds to a F-test prob- 
ability of 4.1 x 10 -4 , approximately a 3.5er significance. Ob- 
servation 1 of 4U 1705-44 shows an decrease of A% 2 = 28.2 
when using a diskline model rather than a Gaussian, which 
corresponds to an F-test probability of 6.6 x 10 -5 , or 4.0cx 
significance. Finally, BeppoSAX observation 1 of GX 17+2 
shows a decrease of A% 2 = 14.7 when using a diskline model 
rather than a Gaussian, which corresponds to an F-test proba- 
bility of 2.3 x 10 -3 , or 3.0cr significance. It is important to also 
take into account the number of observations we searched to 
find the asymmetric lines, as this will reduce the significance. 
In total, we compared diskline and Gaussian fits to a total of 
19 observations from EXOSAT and BeppoSAX. The probabil- 
ities given above should therefore be multiplied by 19, reduc- 
ing the significances to 2.7, 3.2, and 2.0 a for GX 349+2, 4U 
1705-44, and GX 17+2, respectively. 



lEgron et al.l d201 11) discuss how under the criteria of 
Protassov et al. (2002), the F-test ca n be properly appli ed to 
comparing Gaussian and diskline fits. lEgron et ail (1201 1[) also 
apply the posterior predictive p-value test (ppp, iHurkett et alj 
2008) as a method to determine whether a Gaussian or 
diskline provides a better fit. Here, we also apply th is ppp 
test, a nd secondly we perform bootstrap resampling (Efron 
1 19791) of the spectra to further test the significance of the im- 
provement in using the diskline model instead of a Gaussian. 
Essentially, the ppp test involves a Monte Carlo simulation to 
test the likelihood that the diskline model gives the improve- 
ment in x 2 by chance. First, we find the best-fitting Model 
2 using a Gaussian to fit the Fe K line (we use Model 2 as 
it gives better fits than Model 1). Next, we simulate 1000 
sets of spectra (LECS, MECS and PDS spectra for each case) 
with the model parameters randomly drawn using the covari- 
ance matrix of the best-fit and using the detector responses, 
background spectra and exposures from the real data in the 
simulations. These 1000 simulated spectra are then fit with 
the same (Gaussian) model, and also with the diskline model. 
The posterior predictive distribution is defined by the values 
of A% 2 from comparing fits to the simulated spectra with the 
Gaussian and diskline models. The ppp value is then de- 
fined by the number of instances where A\ 2 (simulation) > 
Ay 2 (d ata), see equation 12 and section 3.2 in Hurket tet al.l 
(2008). In Figure [3] we show the posterior predictive distri- 
butions for the three cases (discussed above) where we find 
a significant improvement using diskline based on the F-test 
alone. In all three cases we do not find a single simulation 
with a A% 2 greater than as measured by the data. Given that 
we find zero instances where Ax 2 (simulation) > A% 2 (data) 
we cannot directly calculate a ppp value. However, we can 
infer that the confidence level indicated by the results corre- 
sponds to better than 99.9% level (i.e. less than 1 occurrence 
in 1000 simulations), which indicates an improvement at bet- 
ter than the 3.29cr level. Thus, the simulations strongly sup- 
port the F-test results that a diskline is the preferred model. 

We also perform another test of whether the diskline model 
fits bette r than a Gaus sian by employing bootstrap resampling 
(see, e.g. lEfronTT9 79). This allows an estimation of the distri- 



8 



Cackett et al. 



TABLE 10 

Spectral fits to BeppoSAX observations using a relativistic emission line model. 



Soil I LL' 


Obs N K 

10 22 cm- 2 


Blackbody 

kT (keV) norm T Q (keV) 


Complonized eomponent 
kTfkeV) t 


norm 


r 


3 ower-la\v 


E (keV) 





Diskline 
*m V°1 






EW (eV) 


X 2 (dof) 


GX 349+2 
4U 1705^(4 
GX 17+2 


1 0.51 ±0.02 
1 1.66 ±0.17 
1 2.33 ± 0.09 


1.45±0.02 (6.0±0.2) X 10" 2 0.47 ± 0.01 
1.58±0.O7 (1.8 ±0.3) X 10" 2 0.48 ± 0.03 
1.55 ±0.06 (5.7 ± 0.4) X 10~ 2 0.61 ±0.01 


2.63 ± 0.04 6.7 ±0.2 
2.64±0.12 6.0±0.5 
3.11 ±0.08 5.4±0.3 


1 .18 ± 0.02 
0.47 ± 0.06 
1 .00 ± 0.04 


2.1 ±0.2 
2.9 ±0.1 
3.18 ±0.0 


(5.4+g) x lO" 2 
0.8 ± 0.5 
i 3.2 ±0.7 


697 -0.35 
6 65+ - 32 


-3.0 ± 0.4 

-5 3+ L6 
"-4.7 

< -4.8 


8.3+]° 22+? 
11+J 29 ±6 
18+q 37 ± 8 


(1.2+Oj) x 
(5.6 ± 0.6) X 
(9.0 ± 0.9) X 


O" 2 
10" 3 
lO" 3 


111 ± 13 

120 ± 20 
^7+23 


138.2(103) 
113.7 (87) 
132.0(115) 


Note. — R 


n is given in GM/c 2 


In the diskline model, the outer disk radius is fixed at R ou{ - lOOOGM/e" in 


all cases. 





















1.05 



(c) GX 17+2, Gaussian 



1.05 - 




(e) GX 349+2, Gaussian 




(f) GX 349+2, diskline 




Energy (keV) 



Energy (keV) 



Energy (keV) 



FIG. 2. — Line profiles from BeppoSAX observations of 4U 1705+44 (left, observation 1), GX 17+2 (center, observation 1) and GX 349+2 (right, observation 
1) fit with both a Gaussian (top), and diskline (bottom). In all three cases, the diskline gives a statistically better fit. 



bution in A\ 2 based on the distribution of the observed data. 
Bootstrap resampling consists of resampling the data with re- 
placement. Thus, applying to spectra with N events, we ran- 
domly select N events from the spectrum, with replacement to 
create a new resampled spectrum. In this way some events get 
selected multiple times, where as others will not get selected 
at all. Doing this, we create 1000 resampled spectra from both 
the LECS and MECS detectors. We cannot do this for the 
PDS spectrum, as the pipeline produced spectrum is already 
background subtracted, and thus does not allow for a resam- 
pling of the detector events. However, as the region of interest 
(the Fe K line) is within the MECS detector, this is not an is- 
sue here. Thus, we resample the LECS and MECS spectra, 
but for the PDS we use the original data in the fit each time. 
For each set of resampled spectra, we fit both the diskline and 
Gaussian models (Model 2), and measure the A\ 2 between 
the fits. In Figure [4] we show the distribution of A\ 2 for the 
1000 resamples for each three cases. As should be expected, 
the mean of the distribution is at approximately the A\ 2 of 
the data. Counting the fraction of samplings where Ax 2 > 
(diskline better than a Gaussian) helps understand the signif- 
icance of the diskline being a better fit than the Gaussian. 
For 4U 1705-44 we find 99.4% of resampled spectra have 
A X 2 > 0, for GX 17+2 we find 99.8% and for GX 349+2 we 
find 100%. Again, strongly supporting that a diskline is the 
preferred model over a Gaussian. 

To determine whether the models used are appr opriate, we 
perfo rm a Kolmogorov-Smirnov (KS) test (e.g. Pre ss et al.l 
1992). In the KS tests we determine the cumulative distribu- 



tion function of the BeppoSAX/MECS spectrum as a function 
of energy, and compare this with the cumulative distribution 
function of the model. We do this for both the model with 
a Gaussian and the model with a diskline. We find that for 
all three cases for both the Gaussian and the diskline model 
that the KS probability = 1.0, indicating that both models are 
adequate descriptions of the data. 

It is also important to consider whether we should have ob- 
served more asymmetric line profiles in these archival data. 
Clearly, the spectral resolution of RXTE/PCA is too low to be 
able to show line shapes similar to the Suzaku line profiles. 
For EXOSAT and BeppoSAX, where the spectral resolution 
(FWHM) at 6 keV is 0.6 and 0.5 keV respectively, the answer- 
is not as immediately clear. The asymmetric Ser X-l line as 
seen by Suzaku shows a narrow peak and broad wing which 
is not apparent in the EXOSAT and BeppoSAX data (see top 
left panel of Fig. Q] where the difference between the Suzaku 
and other line profiles can be seen). We have fit the EXOSAT 
and BeppoSAX spectra usin g a diskline with the parameters 
fixed to the Suzaku values (Cackett et al. 2010), except the 
line normalization which was variable in order to match the 
line strength. This intrinsically asymmetric line profile fits 
both spectra reasonably well (x 2 (dof) = 214.9(185) for EX- 
OSAT and x 2 (dof) = 241.1(135) for BeppoSAX with Model 
1), and there are no large residuals in the iron line region, 
though some residuals are present above 7 keV in the Bep- 
poSAX data (we show the residuals in Figure |5}. If the line 
parameters are allowed to be free parameters in the fit the pa- 
rameters are not well constrained for EXOSAT (for instance 
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FIG. 3. — Posterior predictive distributions for testing the likelihood that 
the diskline model gives the observed Ax" by chance for the BeppoSAX 
observations of (a) 4U 1705-44, (b) GX 17+2 and (c) GX 349+2 (ob- 
servation 1 for each object). The distributions are histograms of Ax 2 = 
X 1 (Gaussian) — x' '(diskline) from the 1000 simulations. The red, solid line 
indicates the Ax 2 measured from fitting the data. In none of the cases do 
any simulations show a Ax 2 as large as observed, indicating that the diskline 
model is strongly preferred. 



the inclination tends to 90°). However, for BeppoSAX the \ 2 
improves slightly to x 2 (dof) = 231.5(131), and apart from the 
normalization, the parameters remain close to the Suzaku val- 
ues. For comparison, the best fitting line parameters from fit- 
ting the BeppoSAX data are: £ = 6.97- 004 keV, /3 = -4.3 ±0.5, 
Rin = 6 +1 GM/c 2 , i = 28 ± 1, norm = (9.0 ± 0.7) x 10" 3 . The 
difference between the Suzaku profile of Ser X-l and the EX- 
OSAT and BeppoSAX profiles is the most strikingly of the four 
sources, thus the fact that the Suzaku profile still fits the spec- 
trum reasonably well demonstrates the significant smoothing 
by the lower resolution of those instruments, which may gen- 
erally mask line asymmetry by blurring out the structure. 

Furthermore, we perform additional simulations to see 
whether the A% 2 value that we observe is what would be ex- 
pected if the line profile during the BeppoSAX observations 
is the same as during the Suzaku observations. We there- 
fore take the best-fitting BeppoSAX continuum (Model 2) and 
add a disk line with the same parameters as the best Suzaku 
fits (from iCackett et al.l 120101) . with only the normalization 
changed to match the EW seen in BeppoSAX observation. We 
perform 100 simulations, fitting the simulated spectra with 
both a Gaussian and diskline, and look at the distribution of 
Ax 2 . For GX 349+2, we find that the median A\ 2 = 31, com- 
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FIG. 4. — Distributions of Ax 2 from bootstrap resampling of the Bep- 
poSAX observations of (a) 4U 1705-44, (b) GX 17+2 and (c) GX 349+2 
(again, observation 1 for each object). The distributions are histograms of 
Ax 2 = x 2 (Gaussian) — x 2 (diskline) from the 1000 resamplings. The red, 
solid line indicates the Ax 2 measured from fitting the data. The blue, dashed 
line indicates the mean Ax 2 from the simulations. 

parable to the A% 2 = 22.6 that we observe from the real data. 
For GX 17+2, we find that the median A% 2 = 0.5, much less 
significant than the A\ 2 = 14.7 from the real data. For 4U 
1705-44. we find that the median Ax 2 = 18.4, a little less 
than the Ax 2 = 28.2 from the real data. This demonstrates 
that for GX 349+2 and 4U 1705-44 we see approximately 
the improvement in x 2 th at we would expect based on the 
Suzaku line profiles, where as for GX 17+2 we see an even 
better improvement than expected. 

4.2. Continuum model dependence 

We fitted the data with two different continuum models. 
Model 1 contained a single temperature blackbody, disk- 
blackbody and a power-law (when needed), and Model 2 con- 
tained a single temperature blackbody, thermal Comptoniza- 
tion and a power-law (when needed). We generally find that 
Model 2 gives a better fit to the spectra than Model 1, and 
there is no clear dependence on source state or flux for this. 

We consider how robust the Gaussian line widths are to the 
different continuum models. In Figure |6]we compare the line 
widths (a) obtained when using both continuum models. We 
only compare observations where the two different models 
give equally good fits. We define equally good fits as the 
A\ 2 between the two models being less than 11.8. This is 
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Energy (keV) 

FIG. 5.— Data/model for fits to the (a) EXOSAT and (b) BeppoSAX data 
of Ser X-l. A diskline model was used to fit the Fe line, with all parameters 
ex cept the normaliza tion fixed at those determined from fits to Suzaku data 
in lCacke tt et al. (2010). There are no large residuals in the Fe line region, 
indicating the Suzaku diskline model fits well, though some residuals are ap- 
parent above 7 keV in the BeppoSAX data. Panel (c) is for the same data as 
(b) except that here all the line parameters are allowed to vary. 

equivalent to 3er for the two additional degrees of freedom 
of Model 2 compared to Model 1. The exception to this is 
the Suzaku data where we find that only one observation has 
equally good fits, and so for illustrative purposes we show all 
the Suzaku observations. 

Generally, we find that there is a good agreement between 
the two models when equally good fits are compared. The 
biggest outliers that are not in agreement all have large uncer- 
tainties in the line widths. These are generally EXOSAT and 
BeppoSAX observations - the two missions with the lowest 
effective area, and also significantly lower spectral resolution 
than Suzaku. Given that most observations here have simi- 
lar exposure times, it is clear these observations also have the 
lowest S/N. 

It is interesting to note that the Suzaku and RXTE obser- 
vations almost all give consistent results with differing mod- 
els. While RXTE has the lowest spectral resolution, it has the 
highest effective area, whereas Suzaku has by the far the most 
superior spectral resolution. We can try and understand why 
there is a difference in the Gaussian line width with differ- 
ent models for EXOSAT and BeppoSAX by considering the 
line profiles seen by Suzaku. Firstly, in Figure [7] we show the 
Suzaku line profiles (when fit by a diskline model) for both 
continuum models. Clearly, the line profiles vary very lit- 
tle with the different continuum models, demonstrating that a 
robust line profile can be obtained with high spectral resolu- 
tion and also good S/N. Note that the only Suzaku observation 
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FIG. 6. — A comparison of the Gaussian line widths measured when fitting 
the spectra with two different continuum models (model 1 : diskbb + black- 
body + power-law or model 2: bbody + comptt + power-law). Black circles 
represents data from Suzaku, red stars are from EXOSAT. blue triangles are 
from BeppoSAX, and green squares are from RXTE. 

where there is a large difference in a is the first observation 
of 4U 1705-44, which is the observation with t he lowest S/N 
of all t he Suzaku observations, and one where in Cack ettet al.l 
(2010) we had to fix line parameters because they were poorly 
constrained. 

The Suzaku line profiles generally show a strong and quite 
narrow line peak, with a broader, weaker red wing. A single 
Gaussian does not fit the line profile well and, in fact, two 
Gauss ians (one broad, one narrower) generally provide a bet- 
ter fit (Ca ckett et al.ll2008l) . We also fit the Suzaku data with 
two Gaussians, and give the results in Table QT| (for observa- 
tions where the addition of a second Gaussian improves the 
fit). Figure|8]illustrates how two Gaussians can fit the Suzaku 
line profile ofGX 349+2. 

In the case of RXTE-type spectral resolution, the details of 
the line shape are blurred out by the line-spread function to a 
stage where the lines always appear Gaussian. The intermedi- 
ate spectral resolution of EXOSAT and BeppoSAX, can show a 
hint of a red wing (see panel (c) in Fig. [TJ. However, whether 
a broad or narrow Gaussian fits the line profile best can de- 
pend on the underlying continuum model, though as we note 
above it is the observations with the largest uncertainty in line 
width that are the biggest outliers. 

We therefore conclude that when different continuum mod- 
els fit the data equally well, the line widths determined are 
generally consistent. 

5. DISCUSSION 

We have studied spectra of four neutron star low-mass X- 
ray binaries (Ser X-l, GX 349+2, GX 17+2 and 4U 1705-44) 
from four different missions (EXOSAT, BeppoSAX, RXTE, 
Suzaku) in order to compare the iron line profiles. Our main 
result is that we find iron lines that show very similar line 
profiles between the four different missions examined (Fig. 
Q]). Lines in the CCD spectra of Suzaku are not consistently 
broader than in the gas spectrometer data, demonstrating that 
the broad profiles are not a consequence of pile-up or other 
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FIG. 7. — A comparison of iron emission line profiles when fitting different continuum models to the Suzaku data. Panel (a) shows Serpens X-l, (b) GX 349+2 
(observation 1), (c) GX 17+2 (observation 2), and (d) 4U 1705-44 (observation 2). In all panels, the spectra are from the combined front-illuminated XIS 
detectors. Black is when fitting with a diskbb component, while red is when the comptt component is used instead. Clearly, there is not a strong dependence on 
the choice of continuum model tested here. 



TABLE 1 1 

Gaussian line parameters when fitting two Gaussians to Suzcikti observations with an absorbed diskbb + bbody + power-law 

CONTINUUM MODEL. 



Source 


Obs. 




Gaussian 








G 






x 2 (dof) 






E (keVJ 


er (keV) 






EW (eV) 


E (keV) 


<t (keV) 




EW (eV) 




Serpens X-l 


1 


6.40+0.16 


0.64 ± 0.04 


(3.3 ±0.4) X 10" 


3 


53 ± 10 


6.70 ±0.01 


0.14 ±0.02 


(1.1 ±0.1) X 10 -3 


19 ± 5 


1464.6(1104) 


GX 349+2 


1 


6.40+ - 22 


0.72 ± 0.05 


(1.4±0.1) X 10" 




80 ± 10 


6.79 ± 0.03 


0.21 ±0.03 


(3.1 ±0.5) X 10~ 3 


20 ± 8 


1279.0(1105) 


GX 349+2 


2 


6.40+0.34 


0.73 ± 0.07 


(6.5 ± 1.1) X 10" 


2 


48 ± 13 


6.77 ± 0.02 


0.24 ± 0.02 


(4.0 ±0.5) X 10" 3 


34 ± 10 


1253.5 (1105) 


GX 17+2 


2 


6-40*0-29 


0.50 ± 0.08 


(4.9 ± 1.2) X 10* 


2 


27 ± 11 


6.64 ± 0.02 


0.14 ±0.03 


(2.9 ±0.6) X 10" 3 


17 ± 8 


650.9 (605) 


4U 1705-44 


2 


6.40+ - 12 


0.72 ± 0.04 


(3.3 ±0.3) X 10" 


3 


111 ± 19 


6.70 ± 0.03 


0.12 ±0.03 


(2.5 ±0.6) X 10" 4 


9 ±6 


1338.2 (1120) 



t 1 i 1 r 




Energy (keV) 

FIG. 8. — A illustration that two Gaussians, one broad (green) and one 
narrower (red), can fit the line profile well in the Suzaku data of GX 349+2, 
observation 1. The combination of both Gaussians is shown as a black, solid 
line. 



instrumental effects. Moreover, we found three BeppoSAX 
observations (one of GX 349+2, one of 4U 1705-44 and one 
of GX 17+2) that show evidence for asymmetry, with a rela- 
tivistic disk-line model providing a better fit than a Gaussian. 
However, for data of average quality, the spectral resolution of 
the gas detectors (^0.5 keV for the best case looked at here, 
BeppoSAX/MECS) is not good enough to clearly show asym- 
metric profiles. 

We also found that generally the continuum model choice 
leads to consistent iron line profiles when the models fit 
equally well. However, differences can arise and this is par- 
ticularly a problem with the lower spectral resolution of gas 
spectrometers. However, with CCD spectral resolution and 
high S/N, we demonstrated that the line profiles can be ro- 
bustly determined, regardless of the continuum model. This 
issue is more significant for neutron star LMXBs than black 
hole X-ray binaries due to the typically higher curvature and 
level of continuum degeneracy in the spectra of neutron star 
LMXBs. 
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The extensive simulations of iMiller et al.l (12.010b . the rela- 
tivistic lin e profile detected i n pile-up-free observations in 4U 
1728-34 dEgron et al.ll201 ll) in EPIC/pn "timing" mode, and 
the fact that the Suzaku line pro files do not change regard- 
less of the extraction region used dCackett et al.120101) already 
comprise a large body of evidence suggesting that pile-up and 
instrumental effects are not the source of broad, asymmetric 
lines seen in neutron star LMXBs. Furthermore, our analy- 
sis of archival data presented here also shows that broad lines 
with comparable widths as the Suzaku data are seen in gas 
spectrometer data, where pile-up does not occur. In addition, 
three BeppoSAX observations also show evidence for asym- 
metric line profiles, with a relativistic line model providing 
a better fit than a Gaussian, and with the diskline parameters 
comparable with those observed by Suzaku. 

It is also important to note that inferred upper limits on stel- 
lar radii are consist ent with expectatio ns from dense matter 
equations of state dCackett et al.ll2008l) . Furthermore, mag- 
netic field estimates based on the inner disk radius mea- 
sured from the Fe l i nes in two pulsars, SAXJ1808.4-3658 
dCackett et al.l 120091: iPapitto et alJ I2009T) and IGR J 1 7480- 
2446 in the globu lar cluster Terzan 5 (IMiller et al.l 120111 
Papitt o et alJl201 ll ). are consistent with the magnetic field es- 
timates determined from timing methods. Thus, the measured 
inner disk radii are consistent with expectations. 

The broadest iron lines in neutron star LMXBs, extend 
down to approximately 4.5 keV, i.e. approximately a 2 
keV redshift. If due to Doppler broadening alone, such a 
shift in emission would indicate high velocities of the or- 
der ^0.3c. This, of course, ignores any contribution from 
gravitational redshifts and Comptonization which will also 
broaden the line. Certainly, in highly ionized disks, Comp- 
toni zation will be an imp ortant source of line broadening (see 
e.g. Ross & Fabian 2007), but self-consistent reflection mod- 
eling which include those effects also indicate that the inner 
disk is close t o the neutron star surface where relativistic blur- 
ring i s strong dReis et al.l2009tlCackett et al.l2O10HD'Ai et alJ 
12010) . 

Reflection is further supported by the multip le observations 
of4U 1705-44 considered by lLin et al.l(l201oh . These authors 
show that the blackbody flux and iron line flux are strongly 
correlated in the soft states. This supports the hypothesis that 
the blackbody flux (possibly originating from the boundary 
layer) is the source of irradiating flux leading to reflection, at 
least in the soft states. 

Thus far, studies of iron line variability between the states 
has been inconclusive. In the majority of work that has looked 
at this issue, there is no clear evidence of large changes be- 
tween states dlariaetalj 120091: ICackett et all 120 lot iLin et alJ 
l20"I(it iD'Aiet al.l I2010h . In the hard state of 4U 1705-44 
there has been tentative evidence that the line may be nar- 
rower dLin et al.ll2010t iD'Ai et al.ll2010l) . though this may be 
an ionization effect dReis et alJ2009h . While a thorough study 
of line evolution with state is beyond the scope of this paper, 
it is interesting to point out that we do seem to see variability 
in the line profiles, at least in the case of GX 349+2. 

The line profiles from the three BeppoSAX observations are 
shown in Figure [9] Detailed analysis of these BeppoSAX ob- 
servati ons have previously been presented by Di Salvo et al. 
d2001l) and llaria et al.l (120041) . and both those papers carefully 
looked at variability in the Fe line. The variability we see 
is consistent with their work. iDi Salvo et al.l (120011) found a 
difference in the line EW when comparing spectra from flar- 
ing and non-flaring periods of the first BeppoSAX observation 
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FIG. 9. — A comparison of the line profiles from the three BeppoSAX ob- 
servations of GX 349+2. Observation 1 (black) shows a stronger, broader 
line than observations 2 (red) or 3 (blue). See Table|7jfor spectral parameters 
from these fits. 



- the EW was larger in the non-flaring spectrum when there 
is a ha rd power-law tail extending up to 100 keV. Ilaria et al.l 
(2004) considered all the BeppoSAX observations, creating 
spectra based on the location in the color-color diagram. They 
found a very clear trend showing that the line EW decreases 
with increasing source luminosity. They note that the de- 
crease in EW is due to both an increase in the continuum flux 
and a decrease in the line intensity. Interestingly, the spectra 
where the EW is strongest corresponds to the intervals where 
the hard power-law tail is strongest. Spectra from the flaring 
branch show the weakest line, and that is also where the hard 

p ower-law tail is not d etected. 

ICackett et ail d2010l) and lD'Ai et ail d2010l) both discuss re- 
flection in these neutron star LMXB systems, suggesting that, 
at least in the soft state, the irradiating source of hard X-ray 
flux is the boundary layer. However, in the first BeppoSAX ob- 
servation of GX 349+2, there is also an additional hard pow er- 
law component at higher energies (IDi Salvo et al.ll200Th . If 
this originates from a different region, for example a disk 
corona, then there could potentially be two sources of irra- 
diation with different geometries leading to reflection. 

6. CONCLUSIONS 
To summarize our main findings: 

• The iron line profiles seen in gas-based and CCD-based 
spectra are consistent. This demonstrates that the broad 
profiles are intrinsic to the lines and are not due to in- 
strumental effects (such as pile-up). 

• Several BeppoSAX observations of GX 349+2, 
GX 17+2 and 4U 1705-44 all show evidence for 
asymmetric line profiles. A relativistic diskline model 
fits better than a Gaussian line model as demonstrated 
by the posterior predictive p-value and bootstrapping 
methods. 

• We also found that generally the continuum model 
choice leads to consistent iron line profiles when the 
models fit equally well. However, differences can arise 
and this is particularly a problem with the lower spec- 
tral resolution of gas spectrometers, but that line pro- 
files determined by Suzaku are generally robust to the 
continuum choice. 
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